Nature of excited 0"*" states in ^^'^Gd described by the projected shell model 
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Excited 0"'' states are studied in the framework of the projected shell model, aiming at understand- 
ing the nature of these states in deformed nuclei in general, and the recently observed 13 excited 
0"*" states in ^®*Gd in particular. The model, which contains projected two- and four-quasiparticle 
states as building blocks in the basis, is able to reproduce reasonably well the energies for all the 
observed 0^ states. The obtained B(E2) values however tend to suggest that these 0^ states might 
have a mixed nature of quasiparticle excitations coupled to collective vibrations. 
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Dynamic perturbations of nuclear shapes around the 
equilibrium can give rise to physical states at low to 
moderate excitation energies. Classical examples of such 
motion are (3 and 7 vibrations P, Q > in which nucleons 
undergo vibrations in a collective manner. Traditionally, 
the first excited = 0+ states and the first excited 2+ 
states are interpreted, respectively, as the /? and 7 vibra- 
tional states. While the 2+ collective excitations are bet- 
ter understood theoretically, the nature of the lowest 0+ 
excitation of deformed nuclei still remains under debate 
[E^IEHiQ' The physics of higher 0"*" states is even more 
complex because, on one side, they can predominantly be 
multi-phonon states based on the single-phonons , and 
on the other side, they can be quasiparticle (qp) exci- 
tations in nature. The real situation is, perhaps, that 
the two aspects, collective excitations and qp states, are 
mixed by residual interactions. 

Data on — 0+ states have been relatively sparse. 
In a very recent work by Lesher et. al. , a remarkable 
{p, t) experiment revealed a total of 13 excited O"*" states 
in ^^*Gd, below an excitation energy of approximately 
3.1 MeV. This abundance of 0+ states in a single nucleus 
provides significant new information on the poorly un- 
derstood phenomenon, which has immediately sparked 
off theoretical interest. For this energy range, one may 
think about an explanation through collective modes. In 
fact, Zamfir, Zhang, and Casten 1(1 suggested that many 
of the observed 0+ states may be of two-phonon octupole 
character. Nevertheless, these authors warned also that, 
although the mechanism was excluded in their collective 
models, many of the O"*" states in this excitation energy 
range may be predominantly two-qp in character. 

The purpose of the present paper is to investigate 
whether one can explain the nature of the observed 0'*' 
states in terms of qp excitations. In contrast to Ref. 
[id*!, our calculation here does not emphasize the aspect 
of collective excitation. Although, similar to the work of 
Zamfir, Zhang, and Casten [lOj, we may provide only a 
partial image, our results may shed light on the impor- 
tance of considering the qp aspects in understanding the 
nature of these O"*" states. 

Our study is based on the projected shell model (PSM) 



[lll |. The PSM is the spherical shell model biult on a de- 
formed basis. The PSM calculation usually begins with 
the deformed Nilsson single-particle states at a deforma- 
tion e. Pairing correlations are incorporated into the 
Nilsson states by BCS calculations. The consequences of 
the Nilsson-BCS calculations provide us with a set of qp 
states that define the qp vacuum \(j){e)). One then con- 
structs the shell model bases by building multi-qp states. 
The broken symmetry in these states is recovered by an- 
gular momentum projection [ll| (and particle number 
projection, if necessary) to form a shell model basis in 
the laboratory frame. Finally a two-body shell model 
Hamiltonian is diagonalized in the projected space. 

To determine the deformation at which the shell model 
basis should be built, we first study the bulk properties 
of ^^^Gd including the deformation. In a microscopic cal- 
culation, one searches for energy minima by varying the 
deformation parameter. Here, we calculate the angular- 
momentum-projected energies having the form 



E'{e) = 



{m\Hp'\m) 
{m\p'\m) ' 



(1) 



where is the angular-momentum projection operator 
|l2j | which projects the mean-field vacuum \4>{s)) onto 
states with good angular momentum. As many previous 
PSM calculations for the rare earth nuclei, particles in 
three major shells (N = 4,5,6 for neutrons and N = 
3, 4, 5 for protons) are activated in the present calculation 
for "'^^^Gd. For comparison, unprojected energies 



E{e) = 



{cp{e)\H\cp{e)) 
{<P{eme)) ■ 



(2) 



are also calculated. 

As one can see in Fig. 1, the lowest energy for a 
given angular momentum / is well localized at defor- 
mations varying from e « 0.24 at / — to £ k, 0.28 
at / = 12. The ground-state deformation calculated 
by MoUer et. al. d3] for this nucleus yields e — 0.25. 
These angular-momcntum-projected minima lie at de- 
formations that are slightly larger than the mean-field 
minimum (e « 0.23). Fig. 1 indicates that ^^^Gd is 
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FIG. 1: Angular-momentum-projected energies in ^^^Gd for 
the states with / = to / = 12 as a function of defor- 
mation. The unprojected non-rotating energies (denoted as 
mean-field) are also shown for comparison. 



a stably deformed nucleus against rotation, with a pro- 
nounced energy minimum corresponding to a deformed 
prolate shape. In the following calculations, we thus con- 
struct the shell model basis at the deformation e = 0.26. 
The vacuum state \(f){e = 0.26)) is hereafter written as 

|o). 

Following the spirit of the Tamm-Dancoff method [l^l , 
we build the shell model space by including 0-, 2- and 
4-qp states: 



{|0) 



|o) 



|o) 



^rij ^pj 



|0)} 

(3) 



where is the creation operator for a qp and the index 
n (j>) denotes neutron (proton) Nilsson quantum num- 
bers which run over the low- lying orbitals. Thus, the 
projected multi-qp states are the building blocks of our 
shell model basis: 



(4) 



Here, k labels the basis states and are determined by 
configuration mixing. 

We then diagonalize the Hamiltonian in the projected 
multi-qp states given in Q). In the calculation, we em- 
ploy a quadrupole plus pairing Hamiltonian, with inclu- 
sion of quadrupole-pairing term 



H = Ho 



GMp^P-GgJ^PlPf., (5) 



where Hq is the spherical single-particle Hamilto- 
nian which contains a proper spin-orbit force. The 
quadrupole-quadrupole interaction strength x is deter- 
mined by the self-consistent relation with deformation 
e. The monopole pairing strength Gm is taken to be 



Gm = [20 - 13(iV T Z)/A] I A with "-" sign for neutrons 
and "+" sign for protons. Finally, the quadrupole pair- 
ing strength Gq is assumed to be proportional to Gm, 
the proportionality constant being fixed to 0.20 in the 
present work. These interaction strengths arc the same 
as the values used in the previous PSM calculations for 
the rare earth nuclei 0| . 

In Fig. 2, we show the partial theoretical spectrum in 
^^^Gd up to 4 MeV in energy and 18?i in spin. We empha- 
size that all these states have been obtained by a single 
diagonalization, without any adjustment for individual 
states. Out of these many states, let us concentrate on 
the lowest one at each spin (the yrast band, denoted by 
diamonds), and on all the 0+ states (denoted by filled 
tri-angles). Although it is not a focus of our discussion 
in this paper, we believe that a comparison of the yrast 
band with known data can provide a strict test of the 
model and can provide a useful constraint to the calcu- 
lations of the O'^ states. 

Fig. 3 presents the PSM results for the yrast band in 
^^^Gd (the same values shown in Fig. 2 as diamonds), 
which are compared with the known data (14] . in a plot 
of 7-ray energy versus spin. As can be seen, the data 
are described very well. The calculations predict a sud- 
den drop in the curve at spin / = 16, corresponding to 
a backbending in the moment of inertia. This sudden 
change occurs just at the upper part of the band where 
the current measurement stops. Extension from the cur- 
rent measurement should see this phenomenon and will 
provide a strict test of our model prediction. In this re- 
gard, we note that in the isotonic chain of nuclei, ^^"Dy, 
^^^Er, ^^''Yb, and -'^^^Hf, similar backbending effects have 
been observed and successfully described elsewhere by 
the PSM [Hill. 

We now turn our discussion to the excited 0"*" states. 
Lesher et. al. 13 observed 13 0+ states in the nucleus 
^^^Gd, within an excitation energy range of 1.2 MeV to 
3.1 MeV. In Fig. 4, we plot all the theoretical 0"^ states 
(up to 3.2 MeV in energy; the same values shown in Fig. 
2 as filled triangles) in the order of excitation energies. 
The experimental data [^ are shown in the same plot for 
comparison. The predicted 0+ states are found to be in 
the right energy range, although deviations can be clearly 
seen between theory and experiment. 

What we found impressive is the number of 0"*" states 
predicted by the calculation. The PSM produces a suffi- 
cient number of 0"^ states to be compared with data. As 
described in Eq. @, the total wave function |^m) is a 
linear combination of the (projected) basis states given 
in (jSJ. The basis states in ^ are not arbitrarily selected 
but are taken from all the neutron and proton Nilsson 
orbitals that lie close to the Fermi surface. In ^^^Gd, the 
relevant orbitals are: |[521]i/, |[523]i/, ^[505]i/, §[651]:^, 

f [411]7r, §[413]^, 
In each of these 10 
near-Fermi orbitals, nucleons having opposite signs for 
the K quantum number can couple to a 2-qp state with 
total K = 0. Combination of a pair of 2-qp states can 



^ [642]i^ for neutrons, and ^[411]7r, f [411]7r, 



and 

|[532]7r, and |[523]7r for protons 
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FIG. 2: Theoretical energy spectrum of ^^*Gd calculated up to _B = 4 MeV and I = 18ft. Diamonds are the yrast states and 
filled triangles are excited 0^ states. 
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FIG. 3: Comparison of the PSM calculations and data for 
the yrast states in ^^*Gd in the form of 7-ray energy E^{I) — 
E(I) — E(I — 2) versus spin /. (Data are taken from Ref. 

El). 



further give K = 4-qp states. If one neglects the cou- 
pling of these c[p states to the coUective states, the num- 
ber of low-lying 0+ states in deformed nuclei obtained in 
this way depend solely on the single-particle level density 
and the level distribution near the Fermi surface. Since 
similar conditions can also be found in many other rare 
earth nuclei, we expect such an abundance of 0+ states 
as found in ^^*Gd not to be an isolated case. We pre- 
dict that such an abundance of 0"'" states should also be 
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FIG. 4: Comparison of the calculated 0^ states in ^^*Gd with 
data 13 • 



observed in many other nuclei. 

The large number of 0+ states is difficult to obtain 
within collective models. In Ref. one could obtain 

at most 5 excited 0+ states up to 3.2 MeV in calculations 
with the geometric collective model or the interacting 
boson model (with only s- and d-boson). The reason is 
that within such collective models, the number of degrees 
of freedom of collective motion is limited. Only if one 
considered the odd-parity bosons in an extended boson 
space, could the authors in Ref. E3 S^t more excited 0"'" 
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TABLE I: Predicted 2-qp and 4-qp 0+ states (below 3.25 MeV) in ^^**Gd. 



E (MeV) 


B(E2, 0+ ^ 2+) (W.u.) 


qp-states 


Configurations 


1.004 


1.87 


2-qp 


-f[642]z., 


|[642]z. 




1.321 


0.004 


2-qp 


-^[505]!., 






1.360 


0.014 


2-qp 


-|[411l7r, 


§[411]^ 




1.621 


0.076 


2-qp 


-f [413]7r, 


§[413]^ 




1.636 


0.041 


2-qp 


-f[521]t', 


|[521]t. 




1.716 


0.234 


2-qp 


-|[523]j^, 


|[523]z. 




1.841 


1.102 


2-qp 


-f [532]7r, 


§[532]^ 




2.492 


0.328 


2-qp 


-|[523]7r, 


|[523]7r 




2.631 


0.001 


4-qp 


~^[505]v, ^[505]u, 


-§[411]^ 


,|[411]7r 


2.707 


0.161 


2-qp 


-|[651]i/, 


|[651]i/ 




2.818 





4-qp 


-^[505]u, ^[505]i^, 


-f [413]7i 


,|[413]7r 


2.829 





4-qp 


-|[521]i.,|[521]z/,- 


-§[411]^, 


§[411]^ 


3.028 





4-qp 


-|[521]i.,|[521]!/,- 


-§[413]^, 


f [413]7r 


3.048 





4-qp 


-§[523];/, f[523]!^, ' 


-§[411]^, 


§[411]^ 


3.126 


0.002 


2-qp 


-i[411l7r, 


i[411]^ 




3.168 





4-qp 


-^[505]u, §[521]// 


|[413]7r, 


§[411]^ 


3.192 





4-qp 


-|[521]i/, -|[523]i' 


,|[413]7r, 


§[411]^ 


3.245 





4-qp 


-|[523]i.,|[523]z.,- 


-f[413]^, 


|[413]7r 



states. 

In Table I, we list the 18 calculated 0^ states below an 
excitation energy of 3.25 MeV (slightly higher than the 
highest experimental 0"*" state). Their leading configura- 
tions are also given. Note that each of the configurations 
listed here can not be pure, but is the dominant part in 
the corresponding wave function. One sees that most 2- 
qp states coming from the near-Fermi Nilsson levels have 
lower energies. Due to the varying responses of the resid- 
ual interactions through configuration mixing, one finds 
that some 4-qp states are lower in excitation energy than 
2-qp states. 

Up to now, we have seen that the PSM calculations, 
which explicitly include two- and four-qp states built on 
the basis but no vibrational degrees of freedom, can rea- 
sonably produce a sufficient number of 0"*" states within 
the right excitation energy range. This fact suggests that 
these states are qp states in nature contrary to the work 
of Zamfir, Zhang, and Casten [l^]. However, before we do 
that, we should study the transition properties of these 
states. 

The excited 0+ states can decay to the 2^ state in the 
ground state band through E2 transition. One such a 
transition was measured in Ref. |l7j . Preliminary re- 
sults of Lesher et. al. |3| show 12 such transitions. The 
PSM calculations of the transition from the z-th 0+ state 
to the 2+ state, B(E2, 0+ 2+), are Hsted in Table I. 
Comparing the theoretical B(E2) values with the exper- 
imental data we found that the numbers (in W.u.) 
listed in Table I are one or two orders of magnitude too 
small for many of the transitions. 

We note that the B(E2, 0+ 2+) values (a few 
W.u.) of Lesher et. al. are in average of two or- 



ders of magnitude smaller than the in-band B(E2) val- 
ues of the ground-state band (typically, a few hundred 
W.U.). These small B(E2) values usually indicate that 
the observed 0+ states in ReL ^ should carry signifi- 
cant quasiparticle components. However, our calculated 
B(E2) results are much smaller than the values of Lesher 
et. al., which might indicate an insufficient mixing of 
collectivity to the qp states in our model. This tends to 
suggest that although most, if not all, of these 0+ states 
have significant qp components, they are also mixed with 
correlations induced by the collective motion . These 
correlations can be introduced in the PSM framework by 
inclusion of interactions of higher orders of the multipole 
type [23 , and by addition of the D-pair operators to the 
vacuum state |2l|, which takes both quasiparticle and 
collective degrees of freedom explicitly into account in a 
shell model basis. Generator Coordinate Method, which 
consists of a construction of a linear superposition of dif- 
ferent product wave functions, can also be adopted by 
the PSM. 

In summary, the projected shell model was employed 
to understand the nature of excited 0+ states in deformed 
nuclei. The shell model space consists of projected 2- and 
4-qp states on top of the deformed BCS vacuum state. 
Therefore, the calculation emphasized the quasiparticle 
character of these states, in contrast to previous work 
using collective models. Our results were compared with 
the remarkable example of recently observed 13 excited 
0+ states in ^^^Gd Q- After performing exact angu- 
lar momentum projection and configuration mixing cal- 
culations (by two-body residual interactions) for all the 
possible low-lying qp configurations that can give rise to 
K'^ = 0+ states, we found that the obtained energy lev- 
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els as well as the mimber of the states can reasonably 
explain the data. Preliminary measurements of B(E2) 
values suggest that mixing of these qp states with the 
collective, vibrational motion may not be neglected. 
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